STABLE PHASE RETRIEVAL WITH LOW-REDUNDANCY FRAMES 



BERNHARD G. BODMANN AND NATHANIEL HAMMEN 



Abstract. We investigate the recovery of vectors from magnitudes of frame coefficients 
when the frames have a low redundancy, meaning a small number of frame vectors compared 
to the dimension of the Hilbert space. We first show that for vectors in d dimensions, Ad — A 
suitably chosen frame vectors are sufficient to uniquely determine each signal, up to an 
overall unimodular constant, from the magnitudes of its frame coefficients. Then we discuss 
the effect of noise and show that 8d — 4 frame vectors provide a stable recovery if part of the 
frame coefficients is bounded away from zero. In this regime, perturbing the magnitudes of 
the frame coefficients by noise that is sufficiently small results in a recovery error that is at 
most proportional to the noise level. 

1. Introduction 

Phase retrieval is a topic that is currently extensively researched. Part of the effort is 
directed towards applications in X-ray crystallography, where the Fourier transform dictates 
the form of the measured quantities from which a signal is recovered p3l [8] . It is well known 
that the magnitudes of the Fourier transform need to be complemented by additional infor- 
mation about the signal to make recovery feasible [1], for example the magnitudes of the 
fractional Fourier transform [17] . It is then a challenge whether the additional measurements 
can be realized experimentally. Another main motivation for phase retrieval is quantum 
communication, where quantum states need to be estimated from the relative frequencies 
of outcomes occurring in quantum measurements [15J. In this paper, we investigate the ab- 
stract question of finding a small number of linear measurements such that their magnitudes 
characterize a vector in a finite dimensional Hilbert space, up to an overall unimodular con- 
stant. In addition, we wish to make the recovery procedure resilient against noise affecting 
the magnitude measurements. The central idea is that the redundancy inherent in the frame 
coefficients, resulting from the linear dependencies among the frame vectors, compensates 
the loss of information when passing from frame coefficients to their magnitudes. In fact, 
in frame theory the notion of redundancy is usually understood to be the number of frame 
vectors divided by the dimension of the Hilbert space. In this paper, we choose the frame 
vectors in a specific way to recover the vectors, up to a unimodular constant. We address 
the following main questions: How small can we choose the size of a frame and still charac- 
terize each vector uniquely? What conditions ensure that the vector can be recovered with 
a guaranteed accuracy if the measured magnitudes are affected by noise? 

Several strategies have been applied to the problem of phase retrieval, for example the 
reformulation of recovery in terms of the rank-one hermitian x ® x* |6l HJ [5]. This was 
solved with techniques from compressed sensing by rank minimization in an underdetermined 
system p!Tl [9l [TOl [121 US], or even without rank minimization [21]. However, this technique 
does not specify what redundancy, i.e. the is sufficient for recovery. Other recovery procedures 
embed in even higher dimensional spaces by taking more tensor powers of x with itself [3]. 
Constructions based on expander graphs and the polarization identity give us randomized 
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constructions with an explicit bound on the redundancy [7j. However, this is still far from 
the necessary number of vectors derived from the theory of projective embeddings [16], see 



The recovery procedure outlined here assumes that the signal is realized as a complex 
polynomial of degree d — 1. In the absence of noise, the recovery proceeds in several steps: 

(1) The measured quantities are {\f{ijjj)\'^,\f{ijjj) + uif'^Uj)]"^ : < j < 2n — 2,0 < 
j < 2} where Uj = e^'^*-?/^^""^) and ui = e^'^*'/^. Using the Dirichlet kernel and the 
polarization identity, these are extrapolated to the values and f'{z)f{z) for 
each z in the unit circle. 

(2) From the values of these two functions on the unit circle we determine moments of 
the roots of /. 

(3) The moments determine the polynomial up to an overall unimodular constant. 

We investigate how the presence of noise affects each of these steps, and show that under 
certain conditions, for all sufficiently small e, perturbing the measured quantities up to e still 
gives an approximate recovery with an error of order e. This provable stability only extends 
up to a certain noise level. Numerical experiments show that the domain in which the linear 
error bound holds extends far beyond this level. 

This paper is organized as follows: After fixing the notation, we show in Section [2] that a 
vector, up to a unimodular constant, is determined by a specific choice of 4(i — 4 measure- 
ments. Section |3] establishes criteria for stability of the recovery procedure, complemented 
with the results of a numerical simulation. 



2.1. Definition. The space of complex polynomials of degree at most n is denoted as P„. It 
is equipped with the inner product induced by the Lebesgue measure on the unit circle, so 
p,q & Pn have the inner product 



where the overbar denotes complex conjugation. The space of trigonometric polynomials of 
degree at most n, henceforth called T„, is understood to consist of all linear combinations of 
complex polynomials and of their complex conjugates. 

Thus, Pn is the subspace of analytic functions in T„. On the other hand, the map A : f 
l/P Takes / G P„ to a trigonometric polynomial in T^. The first question we wish to resolve 
is at how many points A{f) needs to be evaluated in order to determine {A/ : |A| = 1}. 
The second problem is that of noisy recovery. If the measured quantities are perturbed, is 
it possible to estimate the set accurately? 

The following theorem shows that to determine a polynomial of degree at most n — 1 up 
to a unimodular constant, is enough to know its magnitude at An — 4 points in the complex 
plane. An essential ingredient is a result of Philippe Jaming's [T7], here the special case for 
polynomials. 

2.2. Lemma. Let d G N. If g & Pd~i then it is determined up to a unimodular constant by 
the values of {gl"^ on two lines L and that intersect in an angle a G M \ vrQ. 

Proof. Without loss of generality we take L = M and Lq = e*"M. Then by the positivity 
of on M, it extends to a polynomial whose roots come in pairs related by complex 
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conjugation. On the other hand, the same apphes to the extension of \g\^ from Lq, to the 
complex plane and reflections of its roots about L^^. Now we have a selection principle based 
on the pairings: We pick as the roots of g the intersection of the sets of roots obtained from 
the two extensions. If g were not determined by this, then it would need to have a root 
in the symmetric difference of the roots from the two extensions. However, the symmetric 
difference is invariant under reflecting first about L and then about L^. This composition of 
the two reflections is an irrational rotation, and so if the symmetric difference is non-empty, 
it gives a dense set of roots in a circle, which means (7 = 0. □ 

2.3. Theorem. Let f{z) = Yl'j=o '^j^'' > ^ '"'^^^ circle as before and Sq, = 4>aiM) U {1} 

with a G R \ vrQ and 4>a{z) = ^^^zj- Then sampling {\f{Zj"^)\'^} on 2d — 1 equidistantly 

spaced points {Zj"^}'^'t^'^ o/ Sq, and {|/(w/)p}^^2^ determines f uniquely, up to an overall 
unimodular factor. 

Proof. We proceed in several steps: 

Step 1. Given / G T^-i, OJ = e^'^*/(^°'~^) and the normalized Dirchlet kernel Dci^i{z) = 
2db"Sj=-d+i^^ then f{z) = Y1]Jq fi.^^)D{zuj''^)- By substitution, for a G C, r > 0, 
/(^) = X^jio^ fi'^ + ruj^)D{{z — a)uj^'-' /r). This means from the values aX2d—l equidistant 
points on a circle we can interpolate any trigonometric polynomial of degree at most d — 1. 
Consequently, the magnitudes of on 2(i — 1 equidistant points on Sq determine 

on the entire circle. Because the circles § and intersect in and u, this also determines 
the magnitudes of / at two of the sample points on §. Once the additional magnitudes 
{|/(u;;)p}^^2^ obtained, the Dirichlet kernel determines the magnitude of |/p on all 
points in S U Sq,. 

Step 2. Using the Cayley map z i-^ and the associated polynomial automorphism 

Wf{z) = {l + z)'-'f[-^^) 

we map both S and to lines L and Lq,, / to a polynomial Wf, and |/p to a trigonometric 
polynomial By conformality, the angle between the lines at is the same as the 

angle between the circles. However, the tangent vector 0^(0) = —{l + u)e^°' has an irrational 
argument, whereas iu does not, so the two circles intersect with an angle that is an irrational 
multiple of vr. By the conformality of the map, the same is true for the intersection of L and 

Lq,. 

Step 3. Next, we use P. Jaming's argument for Wf to show that magnitudes on the two 
lines uniquely determines Wf, up to a unimodular multiphcative constant [T7]. Applying 
the inverse of the map W, the same applies to /. □ 

3. Stable recovery in the presence of noise 

The recovery procedure we outline below heavily relies on the analyticity properties of 
the function space. Although an absolute phase can not be measured, we have access to a 
relative phase such as f'{z)/f{z) for some z & C If these values were known on the entire 
unit circle E> = {z & C : \z\ = 1} then we could recover /, up to an overall multiplicative 
constant from contour integrals and Cauchy formulas as outliend further below. 

We first examine how such integrals are affected by perturbed measurements. Unless noted 
otherwise, for any continuous / : C ^ C, ||/||oo = max^gs 
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3.1. Lemma. Let / : C — ?■ C &e an analytic function, and let pi : S ^ C and • S — ?> M 
with IIpiIIoo < e and ||p2||oo < If there exists a S > such that + P2){z) > S for all 
z G S, then 



ff + Pi f'f 



\f?+P2 l/P 



< 
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Intrinsically the recovery procedure is linked to the moments of the roots of the polynomial 
inside the unit disk. We study how Newton's identities are affected when the moments are 
perturbed. 

3.2. Lemma. Let f he a complex polynomial with Nq roots {zj}^^-^ in the open unit disk. Let 
fi{^) — X1a~o b^^^ ~ rij^il-^ ~ ^j) define the monic factor of f whose roots are precisely the 
roots of f that are inside the unit disk. Given the perturbed moments {jlk}kli of the roots 
such that \ jik — Atfcl < 7 for some < 7 < 1 and fik = Xl^^i f^"^ ^ ^ {!) 2, • • • , Nq} then 
there exists C which only depends on Nq such that {ftk}k=i uniquely determine coefficients 
{bk} with \bk - &fc| < C7 for all k e {1,2, . . . , Nq}. 

Proof. If we knew the values of fik we could recover the actual coefficients using Newton's 
identities, which give the recurrence relation bj^q^^^k = ^\^\=i t^i^No-k+i for all k from 1 to 
Nq. Instead we use our approximations fik to find approximated coefficients hk using the 
recurrence relation 6jVo-fe = ~\Yl!i=ik''if'NQ-k+i with — 1- We inductively show that 
\bk — hk\ is 0(7). For the base case, by assumption 

For the inductive step we note that = | Yl!jli ^\\ — Sj^i — ^0, and if Sk is the set 
of all combinations of k roots of f{z) inside the unit disk, then 

I^JVo-fc 



fj'ibNo - iJiibNol = l/ii - A*i| < 7- 
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Thus, with the inductive assumption that for all j < k, there exists a constant C, such that 



bNo-j - bNo-j 



< Cj^, we have 



bNo-k — bNo-k 



1=1 

k 



< 



No-k+l 
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i^l^No-k+l 



k+l 
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k 
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k 



h ^ \l^lbNo-k+l — l^lbNo-k+l + l^lbNo-k+l — l^lt>No-k+l 
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1 
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No 
— i 

thus Cfe = (l/k) Eli((A^o + l)Cfc„. + suffices. 

Next, we show how the coefficients of the monic polynomial factor containing the roots 
on the inside of the disk can be estimated from perturbed moments. 

3.3. Theorem. Let f{z) = ^^=QO,kZ^, with fixed positive constants m and M' such that 
< m < 1/(^)1 and \f'{z)\ < M' for all z on the unit circle S. Let 

1 

a 



1 + 2(1 + ^) 



and e > with e < am? , pi : S — )■ C and p2 : § — ^ with ||pi||oo < Ib2||oo < e, then 



l^k 

for /c e {1, 2, ... , A^o} observes 

IfJ'k — P-kl 



< 



M' 



7 



(1 — «)m^ \ m 

and if J < I then there exists C which only depends on Nq such that fi{z) = Yl!k=obkZ^ — 
Y\!jli{z — Zj), the monic factor of f whose roots are precisely the roots of f that inside the 
unit disk, has approximate coefficients {6^} with 

e ( ^ M' 



bk - bk\ < C 



(1 — a.)m^ 



1 + 



m 
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Proof. Note that 

e §, (l/P +P2)(^) > m^+p2{z) > - e > (1 - a)m' 
so by the first lemma we have 



f'f + Pi f'f 



I/P + P2 l/P 



< 



1^1 — Q;)m2 



< 



[1 — a)m? 



M' 
m 



If we let A^o be the number of roots of / inside the unit circle, and we let ji^ — Z]j=i ^'ji ^^e 
fcth moment of the inner roots of /, then the residue theorem gives us that for any integer 

k e [0, N] 



1 

2ni 
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Note that Nq — iuLq. Because e < am^ we have that 



1/^0 - /^o| < 



(1 — Q;)m^ 



1 + 



< 
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1 + 
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so rounding flo gives us A^o- Thus by the second lemma, we can recover an approximation 
for fi{z) with approximated coefficients bk such that \bk — bk\ < C7. Now re-expressing 7 in 
terms of e gives the desired result. □ 

3.4. Corollary. If f satisfies the hypotheses of the previous theorem, and in addition \f{z)\ < 
M for z eS, € < ^ for d> N and 



2 1 



(d-l)M+M 



and if g{z) — z^~^f{^), then using the perturbation with pi and p2 as above, we can recover 
an approximation for gi{z) = ^^l^bkZ^ (the monic factor of g{z) whose roots are precisely 
the roots of g{z) inside the unit disk) with approximated coefficients bk such that \bk — bk\ is 
0(e) as e 0. 

Proof. First, we note that ^ = ^on §. Thus, \g{z)\ = \z'^~^f{'z)\ = \ f{z)\ on §. Then because 
we have m < \f{z)\ < M on S, we also have m < \g{z)\ < M on E>. We also know that 
g'{z) = {d-l)lg{z) - z'^-'-y-'il), so that W{z)\ < {d-lMz)\ + \f'{z)\ < {d-l)M + M' 
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on S. Note that on S, = which has a perturbation of P2{z), and 



^{d-l)z\g{z)f-z'-'-'f{z)z'^-^m 

^{d-i)z\m\'-^'mm 

which has a perturbation of (d — l)zp2{z) — z'^piiz). Note that both of these perturbations 
are bounded by de. Thus, g{z) is a complex polynomial that satisfies the requirements of 
the theorem, with replaced by d, M' replaced by (d — 1)M + M', e replaced by rfe, and 
a replaced by /3. Thus by the theorem, using the perturbed functions for /(-z), we can 
recover an approximation for gi{z) with approximated coefficients such that \bk — hk\ is 
0{de) = 0(e). □ 

The objective of the following proposition is to control the error when the recovery of the 
polynomial factors from the inner and the outer roots are combined. 



3.5. Proposition. // / satisfies the hypotheses of the preceding theorem, and in addition 

d 

1 



max;jgs \fiz) \ = M and e < ^ for d > N with 



1 + 2 1 + 



I {d-l)M+M' ^ 



then using the perturbed functions for f{z), we can recover an approximation (up to a 
multiplicative constant) for f{z) = '^k=l,CkZ^ with approximated coefficients c^ such that 
maxjfc |cfe — Cfel is 0(e). 

Proof. Let be the number of roots of f{z) inside the unit disk, and let Ng be the number 
of roots of f{z) outside the unit disk. Note that gi{z) obtained in the previous corollary has 



Zj' 

0. Thus, 



roots (77)^=1 for all roots Zj of f{z) outside the unit disk, with d—l — N additional roots at 



j.w=.''-'-''n(.-i) 

Then if we let fo{z) — z'^''^~^^gi{^), we get 

j=l \ J/ j=l J j=l J 



and so if (2;')^]^ are the roots of f{z) inside the unit disk, by applying the thoerem we get 



3n 



which is a constant multiple of f{z). Thus, we just need an approximation for fo{z)fi{z). In 
terms of the coefficients of gi{z), we have 

No+d-l-N ^ d-l-Ni d-l-Ni 
k=0 k=0 k=0 
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Note that if r is the muhiphcity of the root of gi{z), then for all k from to r — 1, 5^ = 0. 
Thus, the coefficients for the r highest degree terms of fo{z) are equal to 0, and foiz) has 
degree No as it should. Note that we can obtain fo{z) simply by reversing the order of the 
coefficients of gi{z), so our approximated coefficients bk for fo{z) are the same approximated 
coefficients that we obtained in the previous corollary, but in reverse order. Thus \bk — bk\ 
is 0(e) as e ^ 0. In other words, there are numbers Cd,bk which do not depend on e, such 
that \bk &fc| < C*(i,6fcf- Also, from the theorem, we have that the approximated coefficients 
cij for fi{z) = ^^l^o-jZ^ also have error that is 0(e) as e — > 0, so there are numbers C^^aj 
which do not depend on e, such that |a,- — aA < Cda^^- Note that on S 



Ni 



i/.(^)i=ni^-4i<2' 



and 



\fo{z)\ 



No /-, 

^d-l-Ni^N-{d-l) Y\{- - — 



No 

n 



< 2^%ze § 



Thus, the max norms on § of fi{z) and fo{z) are bounded by 2^' and 2^° respectively. 
Because all norms on a finite dimensional space are equivalent, there exists a number K^, 
which non-decreasingly depends only on A^, such that for any complex polynomial h{z) = 
Ylik=o^k'^^ degree less than or equal to iV, we have maxfc<Ar \ck\ < KN\\h{z)\\oo- Note that 



d—l / n 



rd-l-Ni \ f \ 

f{z) = fo{z)Mz) = ( b'i-^-Ni-kz' ]=Y.[J2b^ 

k=0 / \j=0 J n=0 \fc=0 



and thus, for the approximation f{z) — Ylt^ ^kz'' we have that for each k 



Ck — bd-i-Nj-jO-k-j 

j=0 
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Then for each k we have 



Cfe - Cfe 



n n 

= X/ ^d-l-Ni-jak-j - ^2 ^d-l-Ni-jak-j 
j=0 j=0 

n 

— l^d-l-iVi-i^fc-j ~ bd-l-Ni-j^k-j 

n 



— ^ [j^d-l-Ni-jak-j — bd-l-Ni-j^k-j 
j=0 



+ 



bd-l-Ni-jO-k-j ~ bd-l-Ni-jO'k 



(^\bd-i-Ni-j\ i^k-j — ofe-ji + 

j=0 



bd- 



l-Ni-j - bd-l-Ni-j 



I) 



n 

< J2 {{Km 2'''' + Cd,,,_,_,^_^e)Cd,a,_,e + Cd,,,_,_,^_,eK^,2^^) 

n 



j=0 



Thus, since e was bounded above, multiplying two polynomials, whose coefficients have an 
error that is 0(e) gives a polynomial whose coefficients have an error that is 0(e). □ 



Finally, we obtain a finite number of measurements by discretizing and interpolating with 
the Dirichlet kernel. 

3.6. Theorem. Let f{z) = Ylt=l)'^kZ^ be a complex polynomial with degree at most d — 1, 
with fixed positive constants m, M and M' such that m < \f{z)\ < M and \f'{z)\ < M' for 

all z on the unit circle S. Let oo = e^^-i and v = e~ be the {2d — l)th and 3rd roots of 
unity, and let 

1 



1 + 2 1 + 



{d-l)M+M 
m 



'-) 



Let e > with e < j^^^, and assume that we know 2d — 1 values {|/(a;')p + e^^oj^io 



and Qd — 3 values {\f{uj^) + z/-'/'(w')p + with each e^ j < e. Then using only 

these values, we can recover an approximation (up to a multiplicative constant) for f{z) with 
approximated coefficients Ck such that \ck — Ck\ is 0(e). 

Proof. If Dd-i{z) — Sfc=-(d-i) '^'^ normalized Dirichlet kernel of degree d — 1, 

then the set of functions {z i— )■ Dd-iizuj'')}'^^^'^'' is orthogonal with respect to the norm 
on S, and in addition it provides interpolation identity for each trigonometric polynomial 
9{z) = Efc=-(d-i) Cfc^^, 9{z) = Ya=o^^ 9{^^)Dd~i{zu-^). If we let ^o(^) = \f{z)\^ and 
gj{z) = \f{z) + f\z)\'^ for j = 1,2,3, then because each of these is a trigonometric 
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polynomial of degree at most d — 1, we know that we have 



2(d-l) 
1=0 



2(d-l) 



2{d-l) 



ZUJ 



1=0 



1=0 



Let hj{z) = Ym=o ^i,jDd-i(zuj ') be the error obtained when using the Dirichlet kernel 
with the known given values to recover each gj{z). Note that 



\hj{z) 



2(d-l) 

Y ^i,jDd-i{zu;-^) 

1=0 



2{d-l) 



< Y hjDd-i{zuj-')\ < {2d-l)e\Da-i{zuj-')\ = {2d - l)e 



1=0 



Thus we have recovered approximations for the functions gj{z) on S, including go{z) = \ f{z)\'^ 

(2d-l)d ~ d 



with perturbation that is less than {2d — l)e < {2d — 1 ) , ^J^ , , = However, to use the 



previous corollary, we also need an approximation for f'{z)f{z). To obtain this, note that 



3 3 3 

1 ^I7^(^,(;.) + h,{z)) = \ Y'^'SA^) + \Y^'h,{z) 



3=1 



3 



\Y^'\f{z)+i^'f{z)\' + \Y^^h, 



3 



I Y^' {M + ^v'(^)) {f{z) + y^nz)) + ^ X]i7^7i,(^) 

. 3 

3 E +/'(^)/(^) 



1 ^ 



f{z)f{z) + \Yv^h,{z) 



Note that 



\ X]j=i ^^hji^) ^ I Yl^i=i \ hj{z)\ < Thus, we have recovered approximations 



1 v^3 



3 ^3 



for \f{z)\'^ and f'{z)f{z) with error bounded by Then, by using the previous corollary 
with {2d — l)e in place of the e in that corollary, we recover an approximation (up to a 
multiplicative constant) for f{z) with approximated coefficients Ck such that \ck — Ck\ is 
0(e). □ 

3.7. Corollary. Let f be a complex polynomial of degree at most d — 1, let u and v and 

e > satisfy the conditions with m, M and [5 he as in the preceding theorem, then {|/(a;')p + 



^/,o}?io"'^U{|/(a;') + ^V'(^')P + eu}; =0 j=i '^'^th each eij < e determines an approximation 
of f) W to a unimodular constant, with accuracy 0(e). 



,2(rf-l) 3 
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+ e^^olfio 1 Parseval identity, determine 



Proof. The measured quantities ' ^- -^^^'^ 

II /III, up to an error proportional to e. By the point- wise lower bound on |/|, the norm 
is bounded below by ||/||2 > (27rm)^/^, so ||/||2 is also known with accuracy 0(e). Let 

]^Io l&fcP)"^^^ then Cfc = 11/112^^/0 determines an approximation g{z) = Ylk=o^kz'' 

□ 



such that \\f /g — c||oo is 0(e) for some c G C with |c| = 1. 

To illustrate these results we ran simulations of perturbed values to verify the recovery 
procedure. We randomly generated a polynomial, and ran a large number of trials with 
randomized e perturbations on the values needed to apply the theorem. For all sample 
polynomials, we obtained results showing a linear relation between the perturbation e and 
the difference in the coefficients of the recovered polynomial and the original polynomial. 



Even when we perturb the values by 50 times the radius of stability eo 



given by the 



(2d-l)d 

theorem, the linear bound remains intact. This indicates that either the radius of stability 
used in our theorem is not the sharpest value that we can obtain for this, or that unstable 
behavior happens outside of this radius only for pathological examples. 
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Figure 1: Experimentally found recovery error from applying the Newton identities to per- 
turbed measurements. This plot shows that the linear bound for the recovery error in terms 
of the noise level remains intact even for values of e that are much larger than eo. 
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